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ABSTRACT 

We present an investigation into the potential effect of systematics inherent in multi- 
band wide field surveys on the dark energy equation of state determination for two 3D 
weak lensing methods. The weak lensing methods are a geometric shear-ratio method 
and 3D cosmic shear. The analysis here uses an extension of the Fisher matrix frame- 
work to jointly include photometric redshift systematics, shear distortion systematics 
and intrinsic alignments. Using analytic parameterisations of these three primary sys- 
tematic effects allows an isolation of systematic parameters of particular importance. 

We show that assuming systematic parameters are fixed, but possibly biased, re- 
sults in potentially large biases in dark energy parameters. We quantify any potential 
bias by defining a Bias Figure of Merit. By marginalising over extra systematic pa- 
rameters such biases are negated at the expense of an increase in the cosmological 
parameter errors. We show the effect on the dark energy Figure of Merit of marginal- 
ising over each systematic parameter individually. We also show the overall reduction 
in the Figure of Merit due to all three types of systematic effects. 

Based on some assumption of the likely level of systematic errors, we find that the 
largest effect on the Figure of Merit comes from uncertainty in the photometric redshift 
systematic parameters. These can reduce the Figure of Merit by up to a factor of 2 to 4 
in both 3D weak lensing methods, if no informative prior on the systematic parameters 
is applied. Shear distortion systematics have a smaller overall effect. Intrinsic alignment 
effects can reduce the Figure of Merit by up to a further factor of 2. This, however, 
is a worst case scenario. By including prior information on systematic parameters 
the Figure of Merit can be recovered to a large extent, and combined constraints 
from 3D cosmic shear and shear ratio are robust to systematics. We conclude that, 
as a rule of thumb, given a realistic current understanding of intrinsic alignments and 
photometric redshifts, then including all three primary systematic effects reduces the 
Figure of Merit by at most a factor of 2. 

Key words: cosmology: observations - gravitational lensing 



1 INTRODUCTION 

It has recently been shown that the equation of state of 
dark energy could be constrained to a high degree of accu- 
racy using wide and deep imaging surveys (see Albrecht et 
al., 2006; Peacock et al., 2006; for recent and extensive re- 
views). 3D weak lensing has been shown to be a particularly 
powerful way to use the information from such surveys in 
the determination of dark energy parameters (see Munshi 
et al., 2006 for a recent review). 3D weak lensing, in which 
the shear and redshift information of every galaxy is used, 
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has the potential to constrain the dark energy equation of 
state, w^z) = pdc{z)/pdc{z), to Aui(2:) ~ 0.01 using sur- 
veys such as Pan-STARRS (Kaiser et al., 2002) and DUNE 
(Refregier et al., 2006). However, the predictions made thus 
far (Heavens et al., 2006; Taylor et al., 2007) have only in- 
cluded statistical errors and have not included systematic 
effects. Since the scientific goal of many future surveys is 
to constrain Aui(2:) ~ 0.01 such systematic effects have the 
potential to render any cosmological constraints impotent. 

In this paper we will address astrophysical, instrumen- 
tal and theoretical systematic effects relevant to multi-band 
weak lensing surveys in an analytic way. We specifically 
study the systematic effects of photometric redshifts, intrin- 
sic alignments and shear distortion. We consider these three 
primary systematics to have potentially the largest effect on 
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the ability of weak lensing surveys to constrain cosmologi- 
cal parameters. Secondary effects, such as source clustering, 
non-Gaussian effects and theoretical approximations such as 
the Born approximation have been shown to have a smaller 
effect on shear statistics (Shapiro & Cooray, 2006; Scrnboloni 
et al., 2007; Schneider et al., 2002). Note that the effect of 
non-Gaussianity (Semboloni et al., 2007) has only been stud- 
ied via simulations, a full analytic investigation could reveal 
non-Gaussianity to be a larger systematic effect (Takada, 
private communication). Wc will examine the degradation 
that these primary effects may produce in the determina- 
tion of the dark energy equation of state constraints for two 
3D weak lensing methods: 3D cosmic shear (Heavens, 2003; 
Castro et al., 2005; Heavens ct al., 2006) and the shear-ratio 
method (Jain & Taylor, 2003; Taylor et al., 2007). 

The spirit of the approach taken here is to use sim- 
ple, analytic, descriptions of systematic effects. By distill- 
ing complex effects into simple components the change in 
cosmological parameter determination due to any particu- 
lar aspect of a systematic effect can be analysed indepen- 
dently. For example, is the bias or the fraction of outliers in 
photometric redshifts a more important factor? The obvious 
penalty in taking such an approach is that the analytic ap- 
proximations made may not be fully representative of real 
systematic effects. 

The two 3D weak lensing methods are introduced in 
Section 2, however we urge the reader to refer to Heav- 
ens ct al. (2006) and Taylor et al. (2007) for a complete 
and in-depth introduction to the methods. In Section 2 we 
also discuss dark energy parameter prediction. In Section 3 
we introduce the primary systematic effects considered, and 
the parameterisations used. The potential bias in dark en- 
ergy parameters due to each systematic effect is presented 
in Section 4. A marginalisation over systematic parameters 
in presented in Section 5. We conclude and recap with a 
discussion in Section 6. For any technical details concerning 
the 3D weak lensing methods and systematic parameters see 
Appendix A and B. 



2 METHODOLOGY 

To analytically investigate the systematic effect on dark en- 
ergy parameter estimation the Fisher matrix methodology 
will be used throughout this paper. In the case of Gaussian- 
distributed data, which will assumed throughout, the Fisher 
matrix is given for a set of parameters ^ by (Tegmark, Tay- 
lor & Heavens, 1997; Jungman et al., 1996; Fisher, 1935) 

(1) 

where C = S -|- A'^ is the covariance matrix of the method, 
a sum of signal S and noise N terms, and ^ is the mean 
of the signal. Commas denote derivatives with respect to a 
parameter, C,i= dC/d<^i. 

There are a number of methods which employ shear 
and redshift information in order to constrain cosmologi- 
cal parameters. The most widely studied is weak lensing 
tomography, which is an extension of 2D cosmic shear to 
multiple redshift bins, and as such occasionally referred to 
as a 2D-I-1 method. There are various incarnations of the 
tomographic technique. Hu (1999) and Takada & White 



(2004) include all correlations between redshift bins. Takada 
& Jain (2004) introduced bispectrum tomography. Jain & 
Taylor (2003) introduced the concept of taking ratios of to- 
mographic bins which was extended to Cross Correlation 
Cosmography (CCC) by Bernstein & Jain (2004). The ef- 
fect of systematic errors on weak lensing tomography has 
also been extensively studied by for example Ma, Hu & 
Huterer (2006), Huterer et al. (2006), Ishak (2005), King 
& Schneider (2002), Bridle & King (2007), Amara & Re- 
frcgicr (2007). However these papers generally only consider 
each systematic individually here we will combine the effect 
of the three systematics considered. Furthermore this paper 
will focus on two different 3D weak lensing methods. 

2.1 Shear-Ratio 

The shear-ratio method (Jain & Taylor, 2003; Taylor et 
al., 2007) takes the ratio of the average tangential shear 
7* around galaxy clusters in an exhaustive set of various 
redshift bin pairs 

Rii = t- (2) 

By taking this ratio any dependence on the mass, or shape, 
of the lensing cluster drops out resulting in a signal that only 
depends on the geometry of the observer-lens-source sys- 
tem. The cosmological parameters that can be constrained 
are therefore only those that affect the geometry of the 
Universe: Q.m, ^de, wo and Wa- Wo paramctcrisc w{z) = 
Wo -|- Waz/{1 + z) = Wo + Wa{l — d) (Chcvallicr and Po- 
larski, 2001) and we do not assume spatial flatness. The 
redshift range is maximally and exhaustively binned such 
that the bin width at any redshift is equal to the photo- 
metric redshift error at that redshift. The leakage or scatter 
of galaxies between bins due to the photometric redshift 
uncertainty is fully taken into account. The abundance of 
clusters as a function of mass and redshift is modelled us- 
ing the halo model. To constrain cosmological parameters 
(Kitching et al., 2007) the ratio of tangential shears behind 
a cluster is measured directly and fitted by varying the the- 
oretical shear-ratio estimate. The Fisher matrix is therefore 
calculated by varying the mean in equation (1). 

2.2 3D Cosmic Shear 

3D cosmic shear (Heavens, 2003; Heavens et al., 2006) re- 
quires no binning in redshift, describing the entire 3D shear 
field using a 3D spherical harmonic expansion for a small 
angle survey. The transform coefficients for a given set of 
azimuthal £, and radial k [/iMpc"'^], wave numbers are given 
by summing over all galaxies g; 

^{kA) = Y,l'kU{kr')e-''-'\ (3) 

9 

following the conventions of Castro et al. (2005), and we as- 
sume a flat sky approximation. Since the mean signal of the 
coefficients is zero the covariance is varied until it matches 
that of the data (Kitching et al., 2007). The 3D cosmic shear 
covariance depends on the matter power spectrum as well 
as the lensing geometry, so the total parameter set that can 
be constrained is: f2m, f^de, ^b, h, as, wo, Wa, ris and the 
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running of the spectral index a„. Again we do not assume 
spatial flatness, and we have neglected the effect of massive 
neutrinos (Kitching et al, 2008b address the potential of 3D 
cosmic shear in constraining massive neutrinos). We use an 
i'max ~ 5000 and a fcmax = 1-5 Mpc~^ and use the same 
assumptions presented in Heavens et al. (2006). 



2.3 Dark Energy Predictions 

If one has a set of cosmological parameters and a set of 
parameters describing a systematic effect V then the total 
parameter set is given by $ = {0, if>). If extra parameters are 

addcxl to the signal part of a method (either the mean or the 
covariance) the new Fisher matrix becomes a combination 
of the cosmological Fisher matrix F^^ , the derivatives of the 
likelihood with respect to the cosmological parameters and 
the systematic parameter F**' and the systematic parame- 
ters with themselves F'^^ . So that for the total parameter 
set $ the Fisher matrix is defined as 

pee pe^ ■ 



Survey 



DUNE (fiducial) PSl 



(4) 



For all results shown we include a predicted 14-month 
Planck prior (described in Taylor et al., 2007) for which we 
use the parameter set Om, ^de, h, as, ^b, wo, Wa, ris, a„, 
optical depth r and the tensor to scalar ratio r. All dark en- 
ergy constraints quoted arc fully marginalised over all other 
parameters. We assume a ACDM (best fit WMAP3; Spcrgcl 
et al., 2007) fiducial cosmology throughout with = 0.27, 
fide = 0.73, h = 0.71, <T8 = 0.80, Clb = 0.04, wo = -1.0, 
Wa = 0.0, Us = 1.0, a„ = 0.0, r = 0.09 and r = 0.01. 

The performance of a particular survey in terms of its 
ability to measure the dark energy equation of state is com- 
monly quantified in terms of a dark energy 'Figure of Merit' 
(FoM) (Dark Energy Task Force Report; Albrecht et al., 
2006). Paramctcrising the dark energy equation of state's 
redshift behaviour using w{z) = wo -|- Wa(l — a) there exists 
a pivot redshift Zp, at which the constraints using this pa- 
rameterisation minimise, and a corresponding error on wiz) 
at that redshift Aw{zp), this corresponds to rewriting the 
equation of state as w{z) = Wp + Wa{ap — a). The FoM is 
given by the reciprocal of the area of the 1-cr (two parame- 
ter) ellipse at the pivot redshift 



FoM : 



AWaAw{Zp) ' 



(5) 



Note that the Dark Energy Task Force uses the inverse 
area of the 2-a ellipse. The results here will show the ef- 
fect on the FoM for a fiducial survey design, the parameters 
of this survey arc shown in Table 1. We will take this to 
be a 20, 000 square degree survey with a median redshift of 
Zm = 0.9 and a surface number density of no = 35 galaxies 
per square arcminute. We assume a fiducial redshift error of 
a^^z) = 0.025(1 + z) and an intrinsic cUipticity dispersion 
of (Te = 0.247 per shear component (note this differs by a 
factor of V2 from the dispersion on the total shear), this 
is similar to a next generation DUNE-typc experiment (the 
DUNE surface number density though is slightly higher). 
Wo also consider a Pan-STARRS-1 (PSl) type experiment, 
which has substantially different survey parameters that any 
results presented can be tested for consistency over survey 
design, this will be done in Section 6.4. 



Area/sqdeg 


20,000 


20, 000 


^median 


0.90 


0.75 


no / sqarcmin 


35 


5 




0.025 


0.05 




0.25 


0.30 


Table 1. The parameters 


describing the surveys investigated. 


Throughout the fiducial survey will use the 


DUNE survey pa- 


rameters. 
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Awo 0.032 


0.031 


0.021 


Awa 0.177 


0.133 


0.108 


Zp 0.248 


0.295 


0.245 


Aw{zp) 0.015 


0.016 


0.015 


FoM 365 


475 


601 



PSl 


Shear-Ratio 


3D Cosmic Shear 


Combined 


Auio 


0.073 


0.074 


0.029 


AWa 


0.322 


0.288 


0.138 


Zp 


0.310 


0.169 


0.178 


Aw{zp) 


0.028 


0.025 


0.020 


FoM 


110 


140 


363 



Table 2. The baseline constraints for the two 3D weak lensing 

methods for the survey designs considered (described in Section 
2.3). Tlic results shown all include a Planck prior. The FoM is 
given by cciuatioii (5). The combined constraints arc the predicted 
constraints from combining the two methods, discussed in Section 
6.4. 



The predicted baseline constraints for the two surveys, 
for the two 3D weak lensing methods and the combined con- 
straints (see Section 6.4) are shown in Table 2. These con- 
straints are in agreement with other predictions for DUNE 
(for example Amara & Refregier, 2007) and with the Dark 
Energy Task Foce report (Albrecht et al., 2006) for a Stage 
III/IV weak lensing space-based mission for DUNE and 
a Stage II/III ground-based weak lensing survey for Pan- 
STARRS. 



3 PRIMARY SYSTEMATIC EFFECTS 

In this Section we will introduce simple paramotorisations 
that will be used to investigate three systematic effects that 
may have the largest effect on dark energy parameter es- 
timation; photometric redshift estimates, shear distortions 
and intrinsic alignments. 



3.1 Photometric redshift uncertainities 

Methods which can constrain the dark energy equation of 
state necessarily need to include some redshift information, 
since dark energy is an accelerating effect changing the ex- 
pansion history and growth of structure over time (redshift). 
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Note that we will be concerned soley with systematic effects 
inherent in multi-band photometric surveys. Experiments 
such as the Square Kilometer Array (SKA; Blake et al., 
2004) will not have photometric redshift uncertainties, so 
that any systematic effect related to photometric rcdshifts 
could be ignored, however this may be at the expense of 
other potential systematics specific to measuring shear us- 
ing radio data (Chang, Refregier & Helfand, 2004) . 

For wide field and relatively deep surveys, consisting of 
10^ to lO'' galaxies, which can be used to constrain cosmo- 
logical parameters to the percent level, spectroscopic red- 
shifts are currently unfeasible and so photometric redshifts 
must be utilised. There arc many techniques that can be 
used to gain a redshift estimate from photometric data for 
example neural-networks (ANNz; CoUister & Lahav, 2004), 
chi-squared fitting (Hyper-Z; Bolzonella, Miralles & Pello, 
2000) and Bayesian estimation (Benitez, 2000; Edmondson 
et al., 2006; Feldman et al., 2006). However, due to the in- 
herent limitations of the photometric technique, all methods 
result in an error on a given redshift estimation and a scat- 
ter in the relation between the true (spectroscopic) redshift 
and the photometric redshift. The success of a photometric 
redshift estimation procedure is most commonly represented 
as a scattered distribution in the spectroscopic-photometric 
(«s, -?ph) plane. 

In order to retain an independent prediction for the 
systematic effects resulting from photometric redshift esti- 
mation we will present a generic description of the (zs, Zph) 
plane and marginalise over all parameters that are used in 
this description for each 3D weak lensing method. We as- 
sume a Gaussian probability distribution in redshift with 
a fraction of outliers, also with a Gaussian distribution, so 
that the resulting distribution is a sum of two Gaussians 



p{z\zp) 



1 — ./out 



/out 



-(z-Coal^^ph-l-^^bias) /2CT^(2:p) 



-(^-'^cal^ph + ^bia.s) /2['^°(^p)l^ 



27ro-°(«p) 



(6) 



Here we have assumed that the photo-z distribution is cal- 
ibrated, but imperfectly, so that the median spectroscopic 
redshift distribution is biased and inclined so that it lies 
along a line 



■ Ccal^ph "I" Zhias-i 



(7) 



where Zbias is some bias and Ccai is a calibration. A value 
Ccai = 1 would mean that the photometric redshift estima- 
tion is perfectly calibrated to a spectroscopic sample. The 
redshift error a^i^z) is assumed to be unknown and the dis- 
tribution is also assumed to lie between some photometric 
redshift range Grange. We also assume a fraction /out of out- 
lying galaxies in the sample centered at a redshift 2meani 
inclined on a slope described by zZias and c°^x, covering a 
redshift range z°angc- The outlying sample's redshift error is 
also unknown (J°{z). Note that we do not include outliers 
that have low spectroscopic redshifts but a broad range in 
estimated photometric redshift. Our analysis may be pes- 
simistic in this case since by including such a sample some 
redshift biasing effects may cancel-out in this analytic ap- 
proximation. 

The total set of free parameters associated with this 
simple photometric redshift parameterisation is (with fidu- 




Figure 1. A simple description of the (zs, Zph) plane. The dotted 
line is the Zs = Zph around which unbiased, and perfectly cali- 
brated, photometric rcdshifts would lie. The dashed hues arc the 
2s = CcaiZph + ^bias lincs for the main sample (larger block) and 
an outlying sample (smaller block). The fiducial redshift error 
behaviour is az{z) = 0.025(1 -|- z). 



cial values): Zhias = 0.0, Ccai = 1.0, az{z) = 0.025(1 -|- z), 

Zrango = 1.50, /out ~ 0.05, Zmoan = 0.5, ^bias ~ O'lj 

c°^i = 1.0. z°ango = 0.5, a°{z) = 0.05. The fiducial values 
are takcu to be representative of the photometric redshift 
techniques available. The results presented have been tested 
against the fiducial values and there is less than a 2% change 
in the FoM results for a 10% change in the fiducial values. 



3.2 Distortion of the shear 

We analytically investigate the potential systematic effects 
of shear distortions by introducing the following parameter- 
isation of the shear 



7 ■ 



-i'2(p 



7 + 7bii 



(8) 



where we have included an unknown rotation of the shear 
field 0, an uncertainty in the shear measurement A-y and 
a possible bias in the shear 7bias- This parameterisation, 
albeit simplified, should give a good idea of the effects to be 
expected. This is similar to the Shear TEsting Programme 
(STEP; Heymans et al., 2006a) parameterisation where we 
have introduced an extra term due to rotation. The STEP 
parameters are m — 1 — A^, and c = 7bias where <^ = for all 
values of m and c . We take fiducial values of A-y = 1.001 and 
7bias = 0.001, motivated by the STEP results and we use a 
fiducial value = 0.001. This parameterisation represents, 
in an analytical and generic way, the combined effects of 
CCD distortions, instrument effects and inaccuracies in the 
shear measurement procedure due to the shear measurement 
method or image pipeline. We acknowledge that this is a 
simplistic way to include such a variety of complex effects. 
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but in the analytic spirit of this paper this simple model 
should yield a benchmark upon which to gauge how much 
cosmological parameter constraints can be degraded by such 
effects. 

Note that we do not name this 'image distortion' since 
the effects that we aim to parameterise may come from in- 
accuracies in the shear measurement method or instrument 
effects as well as distortions and glitches in the images them- 
selves. 



3.3 Intrinsic alignment efFects 

Intrinsic alignment effects include any non-cosmological, 
contaminating, source of shear. The spurious lensing signal 
from the tidal alignment of close pairs of galaxies (Heavens, 
Refregier & Heymans, 2000; Crittenden et al., 2000; Brown 
et al., 2002; Catelan, ct al., 2001; Heymans & Heavens, 2003; 
King & Schneider, 2003), called the intrinsic- intrinsic (II) 
term, can potentially be removed by ignoring the contribu- 
tion to the shear covariance from close pairs of galaxies to 
any weak lensing statistic. Hirata & Seljak (2004) identified 
a more subtle source of intrinsic shear due to the alignment 
of a background galaxies shear with the tidal field of fore- 
ground galaxies, called the shear-intrinsic (GI) term. Note 
that we will refer to GI and II, the combination of both in- 
trinsic alignment effects will be refered to as lA (Intrinsic 
Alignments). 

We use the fitting formulae given by the numerical sim- 
ulations of Heymans et al. (2006) to the II and GI term. The 
GI fitting formula parameterises the effect using an ampli- 
tude Aqi and a scale dependence 6o (Heymans et al., 2006; 
equation 12) 

(^(z.)e'izi)}e=Eiz.,zi)-^^. (9) 

The free parameters were fitted using n-body simulations. 
E{zs,zi) = D{zi)D{zs — zi)/D{zs) is the lensing efficiency 
of the lens-source pair, where D are angular diameter dis- 
tances. Note that ^gi has units of /iMpc~^ arcmins (see 
Table 4). 

Similarly we use the fitting formula given by the numer- 
ical simulations of Heymans et al. (2006) to the II term. This 
fitting formula parameterises the effect using an amplitude 
An only (Heymans et al., 2006; equation 6) 



{e{za)e*{zb)}e = 



Al 



(10) 



l + {rab/Bu)^' 

Bii = l/i~^Mpc is assumed to be a fixed parameter, where 
Tab is the comoving distance between two galaxies. 

3.4 Overview of the effects of the systematic 
parameters 

For a full description of how the paramcrisation of the sys- 
tematic effects are included in each of the 3D weak lensing 
methods see Appendix A and B. 

To summarise, the possible effect that a systematic can 
have on a method is that it can either introduce extra pa- 
rameters, or add an extra covariance. If the method involved 
is one in which the covariance is varied to match the data 
(as in the 3D cosmic shear case) then an extra covariance 
can also lead to extra parameters. If the method varies the 



mean signal (as in the shear-ratio method) then any extra 
covariance adds noise and no extra parameters. Whether ex- 
tra parameters are added as a result of the covariance is a 
property of the method not the systematic effect. Table 3 
summarises the effect on each 3D weak lensing method of 
including the three primary systematics. It can be seen that 
3D cosmic shear has many more extra parameters, and that 
the shear-ratio method, whilst having fewer extra systematic 
parameters, suffers from more extra covariance terms. 

We note in passing that there exist throe types of sys- 
tematic effect that we define below. We assume covariance 
matrix, C = S + N, is a, sum of signal S and noise N. 

• Type I: systematic ip alters signal (mean of the signal 
or signal covariance; depending on the method) but not the 
noise C = S{'tp) + N. A strong (and correct) prior on 
removes systematic. 

• Type II: systematic ip adds to the covariance but not 
the signal C = S + N + {tp). Strong (and correct) prior 
on '0 docs not remove effect of systematic (errors on param- 
eters arc still increased by ^ even if it is known). ^|) can 
be marginalised over if the method varies the covariance to 
match the data. 

• Type HI: Types I-I-II; which adds both extra parameters 
to the signal and adds a covariance C = S{xIj) + N + (ip). 
A strong (and correct) prior on ip removes the dependence 
of the signal on the systematic parameters but does not 
necessarily remove the extra covariance. 

If we have a Type I systematic effect, i.e. C = S'(V') + N 
which does not add any further covariance then the entire 
systematic effect can be encapsulated by marginalising over 
the total parameter set, which is equivalent to measuring 
(self-calibrating) the systematic parameters from the data 
itself. In this way any degeneracies between the systematic 
parameters and cosmological parameters are taken into ac- 
count. In this case the effect of the systematics could poten- 
tially be alleviated by including extra information (a prior) 
on the systematic parameters. Note that one has to rely on 
any parameterisation used being a good description, the in- 
herent danger is that the form is not accurate enough. 

A Type II systematic effect will add an extra covari- 
ance, C = 5* -|- Ai' -|- C^{ip), which may also contain extra 
systematic parameters, and may also be dependent on the 
cosmological parameter set (or some subset). If the method 
used varies the covariance, as opposed to the mean signal, 
to constrain cosmological parameters then any extra covari- 
ance terms are marginalised over as before. However even if 
extra information is included on the systematic parameters 
there may remain an extra covariance term. For example 
if the covariance C = S + N + AC^{tl') where ^ is a new 
(nuisance) parameter the extra covariance will only be elim- 
inated if ^ = not if the error on A is zero = 0. If 
the dependence of C^{i>) on the cosmological parameter set 
is small this will effectively add an extra noise term to the 
covariance. 

This classification is manifest in Table 3 where we clas- 
sify each of the systematic effects considered in this paper 
for each 3D weak lensing method. Note that the type of sys- 
tematic effect depends on the method not the systematic 
effect itself; the shear-ratio test varies the mean to match to 
the data whereas 3D cosmic shear varies the covariance to 
match the data. Table 3 shows those systematic effects that 
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Shear-Ratio Signal Gov. Type 



pz 


V 


X 


I 


7S 


X 


V 


II 


lA 


X 


V 


II 


3D Cosmic Shear 


Signal 


Gov. 


Type 


pz 


V 


X 


I 


7S 


V 


X 


I 


lA 


X 


V 


II 



Table 3. A summary of the effect of each systematic on the two 
3D weak Icnsing methods. A or X indicates that cither a sys- 
tematic effects the signal or adds an extra covariance term(s). 'pz' 
refers to photometric redshift systematics, '75' to shear distortion 
and 'lA' to intrinsic alignments. The systematics are categorised 
into Types in Section 3.4 



alter the signal (Type I) , those that add an extra covariance 
(Typo II). This highlights the danger of assuming that all 
systematic effects may be removed in future by marginalis- 
ing over extra parameters. 

For a list of all the extra systematic parameters used, 
and their fiducial values see Table 4. The photometric red- 
shift values are interpolated from the literature, for exam- 
ple Abdalla et al. (2007), though no specific reference gives 
definitive values. The shear distortion values are talcen from 
the STEP papers (Heymans et al., 2006a; Massey et al., 
2007). The intrinsic alignment values are taken from an n- 
body simulation of the intrinsic alignment effects, Heymans 
et al. (2006). The GI terms have different values for the 
two methods since the shear-ratio method uses tangential 
shear whereas the 3D cosmic shear method uses the 71 and 
72 components of shear directly. Heymans et al. (2006) give 
different intrinsic alignment systematic parameter values for 
7t, 7x and 7, we take the most likely values of the intrinsic 
alignment parameters. 



4 BIAS IN DARK ENERGY PARAMETERS 

Instead of assuming an extra systematic parameter is mea- 
sured from the data (and marginalising over its value) an 
approach can be adopted in which the extra parameter's ef- 
fect, which will be a bias, on the cosmological parameters is 
estimated whilst the extra parameter's value is assumed to 
be fixed. This bias occurs due to assuming the parameter to 
be fixed, but possibly biased by an unknown amount. The 
cosmological parameter constraints will be smaller than if 
the extra parameter is marginalised over at the expense of 
this bias. 

In Taylor et al. (2007) it was shown how to calculate 
such biases using a Fisher matrix approach. The linear bias 

in a measured parameter S6i due to a bias in a fixed model 
(systematic) parameter Stpj is given by 

Se, = -(F">)T^F^fS^,. (11) 

-F*" is the Fisher matrix of measured (cosmological) param- 
eters and F** is a matrix of derivatives with respect to 
parameters assumed to be fixed and those assumed be be 



Extra Parameter 


Fiducial Value 


Photo-z 




^bias 


0.0 


^calibration 


1.0 


CT;j/(l-h2) 


0.025 


■Strange 


1.5 


Outliers 




/out 


0.05 


bias 


0.1 


calibration 


1.0 


"z 


0.05 


Z° 

'^mean 


0.5 


z° 

range 


0.5 


Shear Distortion 




Aj 


1.001 


<P 


0.001 


Tbias 


0.001 


Shear-Ratio I A (7*) 




j4.Qi/(/i]VIpc ^ arcmin) 




^o/arcmin 


1.32 


An 


0. 15 X lO"'"" 


3D Cosmic Shear lA (7) 




AQi/{hMpc~^ arcmin) 


-1.26 X lO"'' 


60 1 arcmin 


0.90 


An 


0.47 X 10-3 



Table 4. A list of the extra systematic parameters used and the 
Fiducial values adopted. 



measured (see equation 4). We also define where 

m = -C^Stp (12) 

which characterises a systematic parameters biasing effect 
on a cosmological parameter. 

Note that throughout we use to represent the 
marginal error on a parameter and S9 to represent the offset 
in the maximum likelihood value of a parameter. 

To encapsulate the biasing effect in the dark energy 
(w[2p], Wa) parameter space we introduce a Bias Figure of 
Merit (BFoM) which is defined as 

BFoM = \ (13) 

\dw{Zp)dWa\ 

where 5w{zp) and 5wa are the biases in the pivot redshift 
value and the value of Wa due to assuming a systematic 
parameter to be fixed. One would wish to maximise the 
BFoM, its value tending to infinity for zero bias. A desir- 
able bias of less than ~ 0.01 in both w{Zp) and Wa results 
in a BFoM= 10, 000 whereas a poor bias would be of or- 
der BFoM< 100. We stress that this is, in analogy with the 
FoM, a diagnostic tool only so that a conceptual understand- 
ing of the relative effect of the systematic parameters can be 
gained. It does not represent any fundamental aspect of the 
likelihood surface and is of course contingent on both the pa- 
rameterisation of w{z) used and in this case the assumption 
of Gaussianity implicit in the Fisher matrix formalism. Fur- 
thermore it has the potential, as does the FoM to become 
artificially dominated by a good result on one parameter 
masking the poor result of the other. 
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Extra Parameter 


Shear- Ratio 








Sw{zp) 


Swa 


BFoM 


^bias 


-|-U.OO 1 / 


— o.U4if5 


u.ooyy 


^calibration 




— U.Uz / u 


yu.uoo 


CTz 


+0.0530 


-2.4449 


7.7140 


2^range 


-0.0184 


-0.1520 


364.64 


/out 


+0.0047 


—0 3236 


659.12 


zl 
bias 


+0.0037 


-0.1573 


1740.1 


cO 

calibration 


+0.0018 


-0.0560 


9719.8 


^? 


-0.0044 


+0.0539 


4249.7 




+0.0019 


-0.0526 


9997.2 


z° 

"^range 


+0.0018 


-0.0752 


7259.2 


<t> 








Tbias 


— 


— 


— 


-4gi 

^0 

All 


- 
- 


- 

- 


- 

- 


- 


- 


- 


Extra Parameter 


3D Cosmic Shear 








Sw(zp) 


Swa 


BFoM 


^bias 


— U.U04O 


— U.ouoD 


07 01 O 

z ( .y ly 


^calibration 


+U.Uioo 


— U. 


ozy.oZ 


CTz 


in nni ^ 
^u.uu±o 


— U.UOZvJ 


1 9070 
IZU / u 


■Strange 


— U.UUDO 


in 7A7^ 


1 Qc: nr: 

lyo-uo 


/out 


—0.0381 


—0 3588 


73.098 


z° 
bias 


—0 nn9Q 


—0.0349 




C° 

calibration 


n nni ^ 
— u.uuio 


— U.U1C50 


O0C5OZ 




_i_i Q V 1 n— 5 


_i_n nnn4 


1 7 V 1 nS 


Z° 

'^mean 


-0.0022 


-0.0551 


8278.3 


'^range 

1 


-0.0007 


-0.0429 


33846 


-6.1 X 10-s 


-0.0003 


5.8 X 10^ 




+0.0123 


+0.0563 


1440.3 


Tbias 








-4gi 


+0.0084 


+0.0492 


2404.7 


Oo 


-3.7 X 10"'^ 


-2.1 X 10-« 


1.3 X 10^2 


An 


-0.0129 


-0.027 


2830.5 



Table 5. The potential bias on the most likely value of w(zp) and 
Wa due to assuming that an extra systematic parameter is fixed, 
and biased by +0.01. The BFoM is defined in equation (13). A ' — ' 
is shown where the extra parameters do not occur in the signal 
part of the method and as such cannot be marginalised over. 



In Tabic 5 wc show the potential bias in w{zp) and Wa 
due to a +0.01 bias in each extra systematic parameter in- 
dividually, for any given parameter the bias in the others 
is assumed to be zero. All extra covariance and noise terms 
due to each systematic eflFect are included. Note that the 
sign of the bias in w{zp) or Wa is contingent on the sign of 
the bias in the systematic parameter considered (in this case 
+0.01), and the magnitude is dependent on the size of the 
bias considered. If the bias in the systematic parameter was 
larger then the bias in w(zp) or Wa would be proportion- 
ally larger. The values given are meant to be indicative of 
the sensitivity of the 3D weak lensing methods daxk energy 
constraints to each systematic parameter. 

It can be seen that a bias in the photometric redsliifts 
parameter ^bias has the largest effect for both methods. This 
indicates that Zbias would need to be accurate to 1 part in 
10* for the shear ratio method and 1 part in 10^ for the 3D 



cosmic shear method for the most likely value of w{zp) to 
be accurate to ±0.01. 

The shear-ratio method is sensitive to all the photomet- 
ric redshift parameters particularly the bias and calibration 
since these affect the scatter/leakage of galaxies between 
bins at all redshifts (see Appendix A). The method is less 
sensitive to the outlying fraction of galaxies since these have 
a smaller effect on the redshift distribution. 

The 3D cosmic shear method is sensitive to all of the 
photometric redshift parameters, particularly the bias, cal- 
ibration and redshift range. The sensitivity to parameters 
such as the bias and calibration is a result of the parame- 
ters affecting the estimated redshift directly and as such the 
weighting of the shear estimators. The photometric redshift 
extra parameters add uncertainty to the photometric red- 
shift probability distribution. The exact form of the proba- 
bility distribution has an effect on dark energy parameter es- 
timation since, as shown in Castro et al. (2005) and Heavens 
et al. (2006) , the majority of the dark energy signal comes 
from i ~ 1000. Through the Bessel function maxirrmm at 
£ « fermax this corresponds to radial modes of fe « 0.35 
Mpc~^. Photometric redshifts damp the radial modes at in- 
termediate and high k values, at scales of fc > 27r7i/(3000cr2) 
for az « 0.05 this corresponds to fe > 0.03. 

The bias in the shear distortion parameter (approx- 
imately the STEP parameter m) has to be SA-, < 
(O.OlVo.0123) = 0.008 for the bias in 6wp < 0.01. This is 
in approximate agreement with Amara & Rofregior (2007), 
who find the bias in their mo parameter needs to be of or- 
der 0.005, though they focus on investigating the bias in the 
variance a^. 

The cosmological dependence of GI and II terms is small 
(the F^^ term in equation, 11). This low bias for the GI 
and II terms extra parameters suggests that the intrinsic 
alignment terms effectively add extra noise to the 3D cosmic 
shear covariance. In particular the scale dependence of the 
GI term (^o) has a very small potential bias since this only 
changes the overall normalisation of the GI term in a small 
non-linear way (see Appendix B). 

The parameters which have a large bias are also the 
values which should be well measured by the data itself: the 
method is very sensitive to these parameters. In addition 
the parameters which have a very small bias should have a 
small effect on the FoM: the method is not sensitive to the 
parameter and the degeneracy between the extra parameter 
and the dark energy parameters is small. Therefore it is the 
parameters with intermediate values of BFoM that should 
have the largest effect when marginalising over them: the 
method is somewhat sensitive to the parameter of interest 
and there is a large degeneracy between the extra parameter 
and the dark energy parameters. 

We have shown that assuming that certain systematic 
parameters are fixed, but biased, can have a large effect on 
dark energy parameter estimation by biasing the most likely 
values of w{zp) or Wa by a large amount. This suggests that 
marginalising over such parameters must be a more reliable 
way of dealing with such effects. When marginalising over 
systematic parameters error bars on cosmological parame- 
ters of interest will be larger but the most likely value of the 
cosmological parameter of interest will remain unbiased. 

In practice, in order to save computational time, one 
may wish to identify which parameters could cause a bias 
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and then only marginalise over those which appear to be 
troublesome. However in this paper we will continue to 
marginalise over all available parameters so that their po- 
tential effect can be monitored. 

Recently Amara & Refregier (2007) have used the bias 
in cosmological parameters to explore how a survey could 
be designed by fixing the parameter accuracy needed and 
asking what bias could be tolerated that would yield such 
an accuracy. In the case of designing a survey (or photo-z 
method for example) investigating the maximum potential 
bias that can be tolerated [for a given FoM] and then using 
this information as a benchmark upon which to minimise 
the bias in the survey design (or method) is the correct pro- 
cedure. The worst case is that systematic parameters will 
be biased by a large and unknown amount, and one must 
assume this worst case in order to place the most stringent 
constraints on design. 

However in the case of assessing the potential impact of 
systematic effects, on the use of 3D lensing given a survey for 
example, marginalising over parameters which have a large 
biases is the more prudent approach. Instead of assuming 
that parameters are biased by a large unknown amounts, 
the parameter can be marginalised over which takes into ac- 
count the full range of potential values of a parameter, not 
just the largest possible deviation from its fiducial/expected 
value. In this case marginalising over these parameters is the 
best approach; at the expense of a larger error on the cos- 
mological parameters the large bias is negated and the most 
likely value of the cosmological parameter remains intact. 
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Figure 2. The worst-case changes in the figure of merit for the 
fiducial survey with each of the individual photometric redshift 
systematic parameters for the shear-ratio method. The dashed 
line shows the baseline FoM. 



In reality each systematic effect parameter may have prior 
information which can only improve upon the results pre- 
sented in this Section. We investigate the effect of prior in- 
formation on the systematic parameters in Section 6. 



5 MARGINALISING OVER SYSTEMATIC 
PARAMETERS 

Here we investigate the effect of marginalisation over 
systematic nuisance parameters. The systematic parame- 
ters are assumed to be extra parameters that are mea- 
sured/calibrated directly from the data. This will not lead 
to a bias in any cosmological parameters but will increase 
the error through degeneracies with the extra parameters. 

We again use equation (4) to construct a Fisher ma- 
trix containing both the cosmological and extra (nuisance) 
parameters. The resultant predicted marginalised error on 
a given cosmological parameter 9, is given by AO = 

^[_F**]^g^. In this way the new predicted cosmological pa- 
rameter error is marginalised over the predicted systematic 
parameter constraints. This can be compared with the cos- 
mological parameter constraint which does not take into ac- 
count the extra marginalisation over the new (nuisance) pa- 
rameters, M = \J[F'^%g. 

In this Section we will progressively add the primary 
systematic effects to the 3D weak lensing methods in turn 
and asses the impact of the systematic effects on the dark 
energy FoM attainable using each method. We will present 
the effect of each individual systematic parameter, so that 
the source of any reduction in the FoM can be identified, 
aswell as the overall reduction in the FoM due to the sys- 
tematic effects. The results are summarised in Table 6. 

It is important to stress that in this Section we do not 
assume any prior information on the systematic parameters, 
and as such the results presented are a worst-case scenario. 



5.1 Shear- Ratio Method 

As shown in Appendix A and summarised in Table 3 only 
the photometric redshift uncertainties add extra systematic 
parameters to the shear-ratio method, the shear distortion 
and intrinsic alignment effects add extra covariance terms 
which in this case act as extra sources of noise. 

a) Photometric Redshift Systematics 

Figure 2 shows the reduction in the FoM for the various 
photometric redshift parameters for the shear-ratio method. 
The dark energy constraints are most sensitive to marginal- 
ising over a bias in the photometric redshifts. This is due to 
the nature of the effect of Zbias, in which a slight change af- 
fects all redshifts. The dark energy constraints are relatively 
insensitive to the redshift range over which the photometric 
redshifts can be used ( Grange) since this parameter does not 
affect all redshifts and the method is insensitive to the max- 
imum, and minimum, redshifts used. The method is also rel- 
atively insensitive to the outlying population, but note that 
we use a fiducial value of /out = 0.05. The insensitivity to the 
majority of the photometric redshift parameters stems from 
the fact that they affect the detailed form of the photometric 
redshift distribution only which slightly changes the amount 
of scatter/leakage of galaxies between bins. The parameters 
which globally affect all redshifts have a larger effect on the 
FoM. When marginalising over all the photometric parame- 
ters the FoM=163, a factor of 2.2 smaller than the baseline 
FoM. 

b) Photometric Redshift Systematics & Intrinsic 
Alignments 
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By including the intrinsic alignment terms GI and II as 
extra sources of noise the effect on the photometric redshift 
parameters is very similar to the case of including photomet- 
ric redshift systematics alone (see Figure 2). The maximum 
FoM is slightly reduced from the baseline FoM by the in- 
troduction of the GI and II noise terms from 364 to 360 as 
a result of these extra noise terms. The addition of the in- 
trinsic alignment noise terms has a small effect on the dark 
energy FoM. This is due to two reasons, firstly the magni- 
tude of the GI and II terms is at least a factor of 10 times 
smaller than the large scale structure noise terms which af- 
fect this method. Secondly the nature of the GI and II terms 
(Appendix A, equation 21), consisting of four positive noise 
contributions and four negative contributions each from the 
different bin-bin covariant combinations, means that some 
cancellation occurs reducing the effects further. 

c) Photometric Redshift Systematics, Intrinsic 
Alignments & Shear Distortion 

We now include photometric redshift, GI, II and 
shear distortion systematics in the shear-ratio method. The 
change in the FoM with the individual photometric redshift 
parameters, now including the shear distortion, GI and II 
noise terms, is again very similar to when the photometric 
redshift systematics alone are included (see Figure 2). The 
maximum FoM is again reduced from the baseline FoM by 
the introduction of the GI and II noise terms from 364 to 
360, the shear distortion systematic terms have a very small 
effect on the FoM. Marginalising over all the photometric 
redshift parameters and including GI, II and shear distor- 
tion noise terms in the shear-ratio method the FoM= 157 a 
factor 2.2 smaller than the baseline FoM. 

5.2 3D Cosmic Shear 

All three of the primary systematic effects add extra pa- 
rameters to the 3D cosmic shear method, and the GI and 
II intrinsic alignment effects add extra covariance terms. In 
this case, as we progressively add the systematic effects, the 
number of extra parameters will increase. 

a) Photometric Redshift Systematics 

Figure 3 shows the change in the predicted FoM using 
3D cosmic shear due to each of the photometric redshift sys- 
tematic parameters individually. Similar to the shear-ratio 
method the largest effect comes from a bias in the photo- 
metric redshifts, this is due to the sensitivity of the method 
to the photometric redshift distribution, as discussed in 
Section 4. When all photometric redshift parameters are 
marginalised over the resulting FoM is 107, a factor of 4.4 
times smaller than the baseline FoM. 

b) Photometric Redshift Systematics & Intrinsic 
Alignments 

The degradation of the FoM with each photometric red- 
shift parameter is similar with the GI and II covariance 
terms in the 3D cosmic shear method included. However 
the maximum FoM is reduced from the baseline FoM by a 
factor of ~ 3 as a result of the introduction of extra covari- 
ance terms from 475 to 168 (see Figure 4). Marginalising 
over the GI and II extra parameters themselves has a small 
effect. This is due to the relatively poor cosmological depen- 
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Figure 3. The worst-case changes in the figure of merit for 
the fiducial survey with each of the individual photometric red- 
shift systematic parameters for the 3D cosmic shear method. The 
dashed line shows the baseline FoM. 



deuce of the GI and II terms. Hence the GI and II terms 
effectively act as extra sources of noise in the covariance, 
the cosmological dependence of the covariances is small. 

c) Photometric Redshift Systematics, Intrinsic 
Alignments & Shear Distortion 

As shown in Appendix B the systematic parameter 7bias 
only has a second order effect on 3D cosmic shear and since 
the fiducial value of 7bias <S 1 we only marginalise over A^y 
and (p. Figure 4 shows the results of marginalising over each 
of the individual systematic parameters including the in- 
trinsic alignment and shear distortion effects. It can be seen 
that the maximum FoM is again reduced due to the intrin- 
sic alignment terms from 475 to 168. The individual shear 
distortion parameters have a small effect on the FoM as ex- 
pected from the small bias shown in Section 4 implying the 
dark energy constraints have a small sensitivity to these pa- 
rameters. The effects are similar for both parameters since 
for small <f) the parameters' response have approximately the 
same functional form: cos {24>) ~ (l-4(/>2)^2_ Marginal- 
ising over all the systematic parameters results in a FoM= 19 
a factor of 24 times smaller than the baseline FoM. 

We emphasize that the results presented in this Section 
are a worst-case situation, where the systematics are deter- 
mined from the weak lensing data alone. With reasonable 
priors on the systematic parameters (Section 6) the situa- 
tion is markedly improved. 



6 DISCUSSION 

Table 6 shows the FoM, and the pivot redshift error, for each 
of the 3D weak lensing methods after each systematic effect 
is progressively added. Note that we do not display the full 
suite of combinations of systematic effects. 

The baseline constraints from the shear-ratio method 
for the fiducial survey design are shown in Table 2, the 
baseline FoM= 364. Since the intrinsic alignment terms and 
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Figure 4. The worst-case changes in the figure of merit for the 
fiducial survey with each of the individual systematic parameters 
for the 3D cosmic shear method, including GI and II covariances. 
The dashed line shows the FoM before marginalisation over extra 
parameters. 



the overall shear distortion only appear in the noise part 
of the shear-ratio method, the only extra parameters to be 
marginalised over are those from the photometric redshift 
parameterisation. 

For 3D cosmic shear the baseline FoM, using the fidu- 
cial survey design, is FoM= 475. In the case of 3D cosmic 
shear the GI and II terms added extra covariances, and 
provided extra parameters to marginalise over. The shear 
distortion systematic also provides extra parameters to be 
marginalised over. 

The values of bias from Table 5 and the reduction in the 
FoM's shown in Figures 2 to 4 highlight different aspects of 
the relationship between the dark energy cosmological pa- 
rameters and the systematic parameters. For example it can 
be seen from Table 5 that 3D cosmic shear has a smaller 
bias than the shear-ratio method for Zbias. However when 
marginalising over Zbias the reduction in FoM is much larger 
for 3D cosmic shear than the shear ratio method. Figure 
5 explains this apperent discrepancy by showing the con- 
straints from both 3D weak lensing methods in the (zbias, 
Wo) plane. The larger bias in the shear-ratio case is due to 
the smaller degeneracy between Zbias and wo, a shift along 
the degenerate direction of the ellipse projects to a larger 
change in wq for a small change in Zbias- However since the 
projection of the ellipse onto the Zbias axis is small the ef- 
fect of marginalising over Zbias has a small effect on the wo 
constraint. In the 3D cosmic shear case the degeneracy is 
smaller, resulting in a smaller bias, but projection onto the 
2;bias axis is larger resulting in an increase in the wq con- 
straint when marginalising. 

a) Photometric Redshift Systematics 

The effect of photometric redshift and shear distortion 
systematics is approximately the same for both methods. 
The similarity in the overall effects of the systematics on 
the methods, and the very different nature of the methods 
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Table 6. In the worst-case with no informative priors this Table 
shows the FoM for the two methods for the fiducial survey design 
as each systematic effect is progressively added and all systematic 
parameters are marginalised over, 'pz' signifies that photometric 
redshift systematics are taken into account, 'GF and TI' signify 
that intrinsic alignment effects (GI or II) are included, '75' signi- 
fies that shear distortion effects are included. The results shown 
all include a Planck prior. The FoM is given by equation (5). 
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Figure 5. The constraints from 3D cosmic shear (red/dark solid 
line) and the shear-ratio method (green/light solid line) in the 
(^bias! 'i'o) plane marginalising over all other photometric redshift 
systematic parameters, intrinsic alignment and shear distortion 
systematics have not been included. A Planck prior is included. 
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being investigated, means that some general conclusions can 
be made. For the both the shear-ratio and 3D cosmic shear 
methods the effect on the FoM from photometric redshift 
and shear distortion systematics results in a relative reduc- 
tion in the FoM by a factor ~ 2 to 4. 

As shown in Heavens et al. (2006) 3D cosmic shear is 
approximately 10 times less sensitive to a photometric red- 
shift bias than the shear-ratio method. We find again that 
the bias in ■w{zp) in Table 5 due to a bias in Zbias is much 
smaller for 3D cosmic shear than for the shear-ratio method. 
The smaller drop in the FoM when Zbias is marginalised over 
in the shear-ratio method relative to the 3D cosmic shear 
method is due to this larger sensitivity as shown in Figure 
5. This larger sensitivity could be attributed to the binning 
in redshift needed for the shear-ratio method. Any quantity 
(for example a cosmological parameter value) calculated in 
a given bin is calculated assuming that the galaxies are in 
that bin. If there is a bias then the derived quantity will be 
systematically incorrect as galaxies are scattered out of the 
bins having a large effect on the signal. Conversely in 3D 
cosmic shear a bias in redshift merely acts as a slightly dif- 
ferent weighting function in redshift, a slight modification of 
the standard je{r) weighting, so that the shear-ratio method 
is more sensitive to the redshift bias than 3D cosmic shear. 

b) Photometric Redshift Systematics &: Intrinsic 
Alignments 

The effect of intrinsic alignments on the FoM is depen- 
dent on the 3D weak lensing method being used. By includ- 
ing photometric redshift uncertainties and intrinsic align- 
ment effects the FoM can be reduced by up to a factor of 
10, though we stress that this is a worst-case scenario where 
no informative priors have been included. 

The GI and II terms have a small effect on the shear- 
ratio methods dark energy constraints. This is due the GI 
and II contributions to the covariance adding positive and 
negative terms which cancel to some extent (see Appedix A) . 
The GI term is small since we average over a small aperture 
about a cluster, using a larger continous area would increase 
this covariance. The II term is small since the intrinsic- 
intrinsic correlation between any two non-nieghbouring red- 
shift bins is small. 

The reduction in the FoM for 3D cosmic shear is princi- 
pally due to the photometric redshift and intrinsic alignment 
effects. The intrinsic alignment effects reduce the maximum 
FoM by a further factor ~ 2, which is in agreement with 
Bridle & King (2007). Bridle & King (2007) find that for a 
redshift error of <Jz{z)/{l + z) = 0.05 the FoM is reduced by 
a factor of ~ 0.7 by including GI alone, and by ~ 0.6 by in- 
cluding GI and II (from Bridle & King, 2007; Figure 5). This 
is in comparison to what we find, shown in Table 6, that the 
FoM is reduced by a factor of ~ 0.6 by including GI alone 
and by ~ 0.3 by including GI and II. This shows some agree- 
ment between the analyses despite the differences in both the 
3D lensing method investigated and the intrinsic alignment 
parameterisations used. Bridle & King (2007) showed that 
photometric redshifts have to be known to an increased ac- 
curacy when intrinsic alignment effects are included, we find 
a complementary result that when marginalising over pho- 
tometric redshift parameters including intrinsic alignments 
can further reduce the FoM by a factor of ~ 6. We compare 
further with Bridle & King in Section 6.1. It should be noted 
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Figure 6. The constraint from the shear-ratio method in the 
{wo, Wa) plane. The blue (darkest) solid ellipse is the Planck 
constraint, the green (lightest) solid ellipse is the baseline shear- 
ratio constraint alone, and the red (darker gray) solid ellipse is the 
combined constraint without systematics. The outer (black) solid 
line shows the shear-ratio constraint including systematic effects 
and the inner (white) solid line shows this constraint combined 
with the Planck prior. We stress again that this is the worst-case 
degradation. 



that the small effect of marginalising over the GI and II ex- 
tra parameters may be a symptom of the parameterisation 
used. A full investigation of different GI and II paramerisa- 
tions will be the subject of future investigations. 

c) Photometric Redshift Systematics, Intrinsic 
Alignments & Shear Distortion 

It can be seen from Figure 4 and from Table 6 that the 
effect of uncertainty in the shear distortion has a small effect 
however this could be due to the parameterisation used. In 
the case of the shear systematic terms, the values of A^, <p 
and 7bias could be estimated from simulations (as is done 
in STEP; Heymans et al., 2006a; Massey et al., 2007) and 
the shear measurement method could be tuned (i.e. extra 
parameters added to minimise bias) , or constructed ab initio 
as in the case of Miller et al. (2007) and Kitching et al. 
(2008a), to minimise such effects (so that ~ 1, </> ~ 
and 7bias ~ 0). 

For the shear-ratio method, including all three system- 
atic effects and marginalising over all the free systematic 
parameters (in this case only the extra photometric redshift 
parameters), the total reduction in the FoM is a factor of 
2.3 when combined with the Planck prior. Figure 6 shows 
the total systematic effect in the {wo, Wa) plane. 

For 3D cosmic shear the effect of the photometric red- 
shift systematics, shear distortion systematics and the in- 
trinsic alignments in the {wo, Wa) plane is shown in Figure 7. 
By including all three systematic effects, and marginalising 
over all free extra parameters, the FoM becomes FoM= 19, 
a reduction in the FoM by a factor of ^ 20 relative to the 
baseline FoM. 
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Figure 7. The constraint from tlie 3D cosmic siicar mctliod in 
tlie (u!o, Wa) plane. Tlie blue (darkest) solid ellipse is the Planck 
constraint, the green (lightest) solid ellipse is the baseline 3D 
cosmic shear constraint alone, and the red (darker gray) solid 
ellipse is the combined constraint without systematics. The inner 
(white) solid line shows the effect on the combined 3D cosmic 
shear and Planck prior constraint including all systematic effect 
considered in this paper. We stress again that this is the worst- 
case degradation. 

6.1 Including Prior Information on Systematic 
Parameters 

The marginalised results presented thus far have been in 
the self-calibration regime, where the data itself is used to 
measure the extra systematic variables with no priors in- 
cluded. This presents a worst-case scenario in the reduction 
of the FoM. In reality there should exist some information 
on systematic parameters either from different cosmologi- 
cal probes, simulations or data analysis techniques (such 
as photometric redshift code). Table 7 shows the effect on 
the FoM of adding prior information on the systematic pa- 
rameters when photometric redshift, intrinsic alignment and 
shear distortion systematics are included for both 3D weak 
lensing methods. We adopt a Gaussian prior of a given width 
(jp and use the same prior for all systematic parameters. It 
can be seen that in order to recover a substantial proportion 
(> 70%) of the baseline FoM a strong prior needs to be in- 
cluded with Op = 0.001. The relative improvement between 
the two weak lensing methods is very similar. 

For simplicity we use a constant prior for all systematic 
parameters, in reality each systematic parameter will have 
a different prior. For example Zbias can currently be con- 
strained to ~ 2%. Abdalla et al. (2007) have demonstrated 
that using neural net photometric redshift technique (AnnZ) 
the bias can be constrained to ±10~^. This current level of 
constraint is still too large to effectively eliminate the FoM 
degradation, however the results from Abdalla et al. (2007) 
suggest that the IR bands will be able to improve the estima- 
tion of the bias significantly. Results from the Sloan Digital 
Sky Survey (SDSS; Adelman- McCarthy et al., 2005) shown 
that the photometric calibration can be known to 2%. 



I ' I ' I ' I ^ 




Figure 8. The change in the relative FoM with a Gaussian prior 
on Zbiasi assuming other systematic parameters to be fixed for 
the shear-ratio method (green/light, lines) and 3D cosmic shear 
(red/dark, lines). The solid lines are without intrinsic alignment 
systematics and the dashed lines include intrinsic alignments. 

STEP (Heymans et al, 2006a; Massey et al., 2007) 
has shown that shear measurement methods can be well 
calibrated by simulations to within AA-y ~ ±0.01 and 
A7bias ~ ±5 X 10~*. We have shown that marginalising over 
Ay4^ causes a small change in the FoM so that dark energy 
constraints could be robust to distortions due to shear mea- 
surement given some improvement. However if a marginal- 
isation is not done then the shear calibration needs to be 
biased by SA-y < 0.008 (see Section 4). 

Heymans et al. (2006) have probability distributions for 
the intrinsic alignment parameters used here with Ay4ii ~ 
±0.20 X 10~^ and A^o ~ ±0.4 arcminutes. So that the prior 
errors here are overly optimistic for the amplitude of the 
intrinsic alignment terms and too optimistic for the scale 
dependence based on current simulations. However since the 
dependence of the FoM on these parameters is so small this 
should not affect our conclusions. 

6.2 Prior on Photometric Redshifts 

It has been shown that the largest effect from a free sys- 
tematic parameter on the FoM comes from the bias in the 
photometric redshifts for both 3D weak lensing methods. We 
should therefore expect that even adding prior information 
only on this parameter would improve the FoM. Figure 8 
shows the improvement in the relative FoM if a Gaussian 
prior is added to the Zbias constraint, marginalising over all 
other parameters which have no added prior. It can be seen 
that ~ 80% of the FoM can be recovered using a Gaussian 
prior with an error of (7p(2:bias) ~ 5 x 10"* for the shear- 
ratio method and (Tp(2:bias) ~ 2 x 10"'^ for 3D cosmic shear 
when intrinsic alignment effects are included (dashed lines). 

To begin to improve upon the FoM for which no prior 
has been added the shear-ratio method requires a prior with 
an error that is approximately 10 times smaller than for the 
3D cosmic shear method. This factor of 10 between the two 
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Shear- Ratio FoM/FoMmax 


3D Cosmic Shear FoM/FoMmax 


Combined FoM/FoMmax 


Self Calibration 


0.43 


0.11 


0.83 


o-p = 0.01 


0.44 


0.37 


0.89 


up = 0.001 


0.60 


0.90 


0.92 



Table 7. The effect on the predicted FoM of adding a Gaussian prior information on the systematic parameters for the two 3D weak 
lensing methods. The self-calibration has no prior information and is a worst case, crp represents a constant Gaussian prior of a given 
width on all the systematic parameters. FoMmax is the maximum FoM given that all three systematic effects are included; for the 
shear-ratio method this is FoMmax = 360 for 3D cosmic shear FoMmax = 168 and for the Combined constraint FoMmax = 594. These 
are reduced from the absolute maximum, baseline FoM, by the inclusion of the intrinsic alignment systematics. 



methods reflects the results of Section 4 and Heavens et al. 
(2006), the 3D cosmic shear method is less sensitive to this 
parameter so that in the self-calibration regime, with a poor 
prior, the reduction in the FoM is larger. Since the constraint 
on Zbias is already smaller for the shear-ratio method, as can 
be seen from Figure 5, the extra prior error needs to be even 
smaller to begin to improve upon the constraint. 

The solid lines in Figure 8 show the improvement in 
the FoM as the prior on Zbias is improved where intrinsic 
alignment effects have not been included. For the shear-ratio 
method this has little effect, since the method is relatively 
insensitive to intrinsic alignment effects as discussed in Sec- 
tion 6. For 3D cosmic shear the requirement on the accu- 
racy of Zbias to recover ~ 80% of the FoM is relaxed to 
c'"p(2bias) ~ 6 X 10~^, a factor of ^ 3 times larger than when 
intrinsic alignments are included. This is in agreement with 
the results of Bridle & King (2007) who find that the average 
photometric redshift error needs to be ~ 3 to 4 times smaller 
to recover ~ 80% of the FoM when intrinsic alignements are 
included. 

Comparing the two 3D shear methods, we find that the 
3D cosmic shear has slightly better ideal FoM, but this de- 
grades more if the systematics need to be estimated from 
the data themselves. The better the prior is on the system- 
atics, the better the 3D cosmic shear method will fare, but 
with reasonable priors which should be achievable with cur- 
rent techniques, the shear-ratio method and 3D cosmic shear 
should attain comparable accuracy. 



6.3 Spectroscopic Calibration of Photometric 
Redshifts 

Spectroscopic redshifts could be used to calibrate the red- 
shift bias. Assuming Poisson statistics the number of spec- 
troscopic redshifts required can be written as 



Combined : pz + lA + 7g 



Cr(2bias) 



(14) 



Using the result from Taylor et al. (2007) this can also be 
related to the bias parameter, defined in equation (11), by 



Cbiasg-z(z) 



Sw{Zp) 



(15) 



For the shear-ratio method Figure 8 shows that a prior er- 



ror of (7p(2:bia 



5 X 10 is needed to eliminate the ef- 



fect of marginalising over Zbias- Assuming an average red- 
shift error of (Tz[z) ~ 0.025 and using equation (14), this 
implies A'spcc ~ 3 x 10"^. If Zbias is assumed to be fixed 
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Figure 9. The change in the figure of merit with each of the 
individual systematic parameters combining the shear-ratio and 
3D cosmic shear methods, including photometric redshift, shear 
distortion, GI and II systematic effects. The dashed line shows 
the baseline combined FoM. 



but biased then from Table 5 we have Cbias ~ 55 for the 
shear-ratio method. If the bias in w{zp) is required to be 
less than 5w{zp) ~ 0.01 then equation (15) implies that 
A'apec ~ 2 X 10'*. These predicted numbers of spectroscopic 
redshifts is in agreement with predictions for weak lensing 
tomography, for example Abdalla et al. (2007). 

For 3D cosmic shear Figure 8 shows that a prior error 
of (Tp(2bias) ~ 2 X 10"'^ is required to recover the original 
FoM, using equation (14) this implies the number of spectro- 
scopic redshifts needed to be N^pec ~ 6 x lO'^ . Assuming Zbias 
is fixed it can be deduced from Table 5 that Cbias ~ 5.48 
which implies, using equation (15) that Nupcc ~ 2 x 10^. 
This relatively small calibrating sample is in approximate 
agreement with the results of Heavens et al. (2006). This is 
as a result of the small degeneracy between Zbias and wo as 
discussed in Section 6. 



6.4 Combined constraints &: Bottom-Line 
Predictions 

By combining the constraints from the shear-ratio and 3D 
cosmic shear method the relative decrease in the FoM may 
be smaller due to parameter degeneracies between the sys- 
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3D CS SR COMB 

Figure 10. The change in the FoM for a DUNE-type survey using 
3D cosmic shear (3D CS), the shear-ratio method (SR) and the 
combined FoM (COMB) for each of the optimistic, intermediate 
and pessimistic scenarios. Here the absolute FoM is shown for the 
fiducial survey considered. 
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3D CS SR COMB 

Figure 11. The change in the FoM for a PSl-type survey using 
3D cosmic shear (3D CS), the shear-ratio method (SR) and the 
combined FoM (COMB) for each of the optimistic, intermediate 
and pessimistic scenarios. Here the absolute FoM is shown for the 
fiducial survey considered. 



tematic parameters being lifted, and the baseline FoM will 
be larger as the dark energy constraints are combined. 

The combination presented here does not take into ac- 
count the full covariance between the two methods but 
should be valid to a first approximation since the shear- 
ratio uses clusters (small scale features) and the majority of 
the dark energy signal for the 3D cosmic shear method is 
from approximately sub-degree scales (the maximum signal 
is at f ~ 1000; Heavens et al., 2006). Also the full correla- 
tion should only appear between the noise of the shear-ratio 
method and the signal of 3D cosmic shear both of which 
depend on the matter power spectrum. 

The baseline FoM for the fiducial survey when the meth- 
ods are combined, with a Planck prior included, is FoM= 601 
with pivot redshift error of Aw(2;p) = 0.015. Figure 9 shows 
the effect of the combined FoM on the systematic parame- 
ters. It can immediately be seen that the combined FoM is 
substantially larger than either method alone, and that the 
degradation of the FoM with the systematic parameters is 
much smaller. In particular the methods' different parameter 
degeneracies are very complementary for ^bias and Ccai . Since 
the shear-ratio method provides no constraint on the shear 
distortion parameters these have a relatively large effect on 
the combined FoM via the 3D cosmic shear method's depen- 
dence. When marginalising over all systematic parameters 
the FoM= 494 with a pivot redshift error of A.w{zj,) — 0.018 
with a Planck prior included. 

This is an encouraging result, despite the naive addition 
of the two methods, it appears that the 3D cosmic shear 
and shear-ratio methods are very complementary in terms of 
both cosmological and systematic parameter degeneracies. 

In Table 8 and Figures 10 and 11 we present the effect 
on the relative and absolute FoM given three different fu- 
ture scenarios, now showing results for both a Pan-STARRS 
survey and the DUNE-type fiducial survey. 

Pessimistically one could assume that the II and GI ef- 



fects cannot be removed from data and that there are no 
prior constraints available for photometric redshift parame- 
ters. This worst-case seems unlikely since already using cur- 
rent photometric redshift codes one could provide priors for 
the bias and calibration of a photometric redshift sample 
(e.g Abdalla et al, 2007) and there are planned and ongo- 
ing spectroscopic instruments that could provide adequate 
calibration. This result is a worst-case scenario, since it has 
been shown (for example Heymans & Heavens, 2003; Hey- 
mans et al., 2004; King & Schneider, 2003; King, 2005) that 
intrinsic alignments can be removed to some extent from 
cosmic shear data. The worst-case FoM~ 20 for the 3D cos- 
mic shear is approximately the same for both surveys so 
that the relative change for DUNE is larger, in this regime 
so much information is lost through systematic effects that 
the differences in survey design have a sub-dominant effect. 

The intermediate stage takes into account current be- 
lief in the likely removal of the II intrinsic alignment term 
and includes a currently realistic prior on the photometric 
redshift systematic parameters; this is currently the most re- 
alistic scenario. There is also broad agreement between the 
relative FoM changes given the different survey designs. We 
conclude that such a scenario could result in a factor of ~ 2 
reduction in the statistical FoM for both 3D weak lensing 
methods. 

Going beyond current photometric redshift constraints 
and understanding of intrinsic alignments we present an op- 
timistic conclusion based on the assumptions that both GI 
and II effect could be removed and the photometric redshift 
systematic parameters could be calibrated to 1 x 10~^ . Again 
there is broad agreement over survey design. In this case the 
3D cosmic shear constraints are degraded by ~ 1.2 and the 
shear-ratio constraints show a FoM reduction by a factor of 
~ 1.7 for the fiducial survey. By combining the 3D cosmic 
shear and shear-ratio constraints parameter degeneracies are 
lifted to such an extent that even in the pessimistic case the 
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reduction in the FoM due to all three primary systematics 
effects is only a factor of 1.2. 

Given the agreement between two different 3D weaJj 
lensing methods and over two different survey designs we 
conclude that one may expect at most a factor of 2 reduc- 
tion in the FoM using 3D weak lensing as a result of the three 
primary systematic effect considered in this paper. We ex- 
pect that this factor would be substantially less in actuality 
given that shear-galaxy correlations could be used and the 
weak lensing methods can be combined as shown. 



7 CONCLUSION 

We have shown using simple analytic approximations that 
systematic effects can have a substantial impact on the abil- 
ity of 3D weak lensing methods, the shear-ratio method and 
the 3D cosmic shear method, to constrain the dark energy 
equation of state. We used the Figure of Merit (FoM) to 
gauge the ability of a next generation experiment to con- 
strain the dark energy equation of state. The systematic 
effects we considered are those associated with photomet- 
ric redshifts, an overall distortion in the image plane and 
intrinsic alignments (both the GI and 11 terms). 

The dark energy FoM can be degraded a factor of 2 
due to the photometric rcdshift systematics alone, for both 
3D weak lensing methods. Shear distortion systematics have 
a small effect for both methods. Intrinsic alignment effects 
alone can degrade the FoM by a further factor of 2 for the 3D 
cosmic shear method, but have a small effect for the shear- 
ratio method. This difference is due to the way in which the 
methods use shear information. 

When an extra systematic parameter is encountered 
it can either be marginalised over using the available 
data (equivalent to self-calibration), thereby increasing the 
marginal error on any cosmological parameter of interest, or 
it can be assumed to be fixed. If a parameter is fixed then any 
deviation away from the assumed value will bias the most 
likely value of any measured cosmological parameter. This 
bias is a function of both the cosmological parameter error 
and the sensitivity of a method to any systematic parameter. 
From this analysis it hcis been shown that assuming some 
parameters to be fixed can lead to large biases in w{zp) and 
Wa, this is complimentary to the analysis done by Amara 
& Refregier (2007) who investiagted the bias in cosmologi- 
cal parameters due to shear measurement systematics using 
weak lensing tomography. 

The methods are remarkably insensitive to many pa- 
rameters, in particular most of the photometric redshift pa- 
rameters, including the fraction of outliers, and an overall 
distortion of the shear field. By adopting a parameterisar 
tion of the photometric-spectroscopic plane we have shown 
that a bias in any photometric redshift redshift technique 
needs to be known to within ±10^'* for the dark energy FoM 
to remain unaffected. We have shown that to calibrate the 
photometric redshift parameters approximately 10^ spectro- 
scopic redshifts are required for the shear-ratio method and 
approximately 10'^ to 10^ for the 3D cosmic shear method. 
This difference can be attributed to the binning in redshift 
required by the shear-ratio method, which may also explain 
the agreement between the predicted spectroscopic require- 
ments of the shear-ratio and shear tomography methods. 



The intrinsic alignment terms were modelled using the 
flcymans et al. (2006) analytic approximations. The Gl and 
11 terms had a small effect on the FoM from the shear-ratio 
method, we found a drop in the FoM of approximately 10% 
when the extra covariances were included. For 3D cosmic 
shear the FoM is reduced by approximately 50%, but there 
is a very small sensitivity to the extra intrinsic alignment 
systematic parameters. The 3D cosmic shear result is in 
agreement to what has been found using shear tomogra- 
phy in Bridle & King (2007). The caveat to these compar- 
isons is that we use a fully 3D cosmic shear method, with 
no binning, and only investigate the Heymans et al. (2006) 
parametcrisation; Bridle & King (2007) consider a variety 
of parametcrisations and investigate a binning tomographic 
method. 

We have shown that a good prior on systematic param- 
eters can improve the FoM. The relative reduction in the 
FoM can be limited to < 30% if the prior on all systematic 
parameters has a Gaussian error of ap = 0.001. In partic- 
ular we have shown that a prior on «bias can improve the 
FoM. Good priors on ^bias, Ccaiibration would be particularly 
helpful in limiting the effect on the FoM. 

By combining the 3D weak lensing methods the FoM 
can be increased by a up to a factor of 6 relative to the meth- 
ods individually. Furthermore the photometric redshift sys- 
tematic parameter degeneracies axe complementary leading 
to less systematic degradation in the combined constraints. 

The bottom line is that the most important systemat- 
ics to control are those concerning the photometric redshift 
distribution. If these can be controlled to 1% accuracy, then 
the FoM for proposed future surveys such as DUNE and 
Pan-STARRS may be reduced by at most a factor of order 
two. In order to reduce these systematics to a negligible level 
needs 0.001 accuracy in median redshifts, requiring 0{W^) 
redshifts. We make a number of recommedations and obser- 
vations with which to guide future systematic investigations 

• Photometric redshift systematics play a dominant role 
in the systematics that affect 3D weak lensing. The indi- 
vidual systematic parameter which can have the largest ef- 
fect on the FoM is the bias in photometric redshifts. How- 
ever approximately 0{W^) spectroscopic redshifts should be 
enough to calibrate photometric redshifts to the required ac- 
curacy. These would need to be representative of the photo- 
metric galaxies and complete. 

• If shear calibration bias is assumed to be fixed then 
an uncertainty in the bias of ~ 0.008 can bias dark energy 
parameters by > 0.01. However marginalising over shear bias 
has a smaller effect on the FoM. 

• Intrinsic alignments play a major role in 3D weak lens- 
ing systematic effects, and can reduce the maximum achiev- 
able FoM by up to 4. The broad agreement between the 
parameterisations investigated here and in Bridle & King 
(2007) suggest that the general trends are robust. 

Despite degrading systematic effects 3D weak lensing 
retains the potential to be the most powerful cosmological 
probe of dark energy. Weak lensing is entering a formative 
period in its development, given that the statistical ability 
of the method to constrain cosmology is accepted attention 
must now be focussed on understanding and reducing sys- 
tematic effects. 
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Survey 




lA 


Phot;0~z Prior 




'\Ti PC! FnAA /ThViA/T 
OU KjiD rOLvL/ iroivimax 


v^omu roiviy jroivimax 


DUNE (fiducial) 


Optimistic 


None 


<Tp = 0.001 


0.58 


0.85 


0.86 


PSl 


Optimistic 


None 


ap = 0.001 


0.76 


0.65 


0.94 


DUNE (fiducial) 


Intermediate 


GI only 


ap = 0.01 


0.45 


0.59 


0.84 


PSl 


Intermediate 


GI only 


o-p = 0.01 


0.56 


0.54 


0.92 


DUNE (fiducial) 


Pessimistic 


II and GI 


None 


0.43 


0.04 


0.82 


PSl 


Pessimistic 


II and GI 


None 


0.44 


0.11 


0.83 



Table 8. Here we present three possible scenarios representing optimistic, intermediate and pessimistic points of view on the future 
development of weak lensing systematics using 3D cosmic shear (3D CS) , the shear-ratio method (SR) and the combining the methods 
(COMB) The pessimistic scenario is unlikely given that II effects can be removed from data and the current photometric redshifts have 
at least some prior on their bias and calibration. Note here FoMmax is the baseline FoM given in Table 2, not the reduced maximum 
FoM used in Table 7 which already included in intrinsic alignment systematics. 
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APPENDIX A 

In this Appendix we will present the technical details of how 
the primary systematic effects are included in the shear-ratio 
method. For full details on the theory and implementation 
of the method see Taylor et al. (2007) and Kitching et al. 
(2007). 



So any extra parameters used to describe the photo- 
metric redshift probability distribution enter the signal of 
shear-ratio method via the determination of the mean tan- 
gential shear in a given redshift bin. 



Shear-ratio photometric redshift systematics 

The shear-ratio method takes the ratio of average tangential 
shear in different redshift bins behind a cluster. The binning 
is necessary since if it is assumed that observed ollipticitics 
have zero mean and that the distribution of intrinsic elliptic- 
ities is Gaussian the resulting distribution of the ratio of el- 
lipticities has a Cauchy/Lorentzian distribution (and hence 
infinite variance). Around a cluster the mean ellipticity is 
non-zero and so Gaussian errors can be assumed, and this is 
certainly so when binning around a cluster and in redshift, 
so that the number of galaxies contributing to the mean is 
increased. 

Photometric redshift systematics affects the shear-ratio 
signal via the uncertainty in a galaxies position meaning it 
could be scattered out of (or in to) a bin i.e. the galaxies 
used to infer the tangential shear in any given bin could in 

actuality be in another bin. The average tangential shear in 
the i'' 

by 



redshift bin behind a cluster at redshift zi is given 



(■yt,i) = 7t,o 



i: 

L 



Sk[r{z) - r(2i)] 
dz — „ , , — —n{z)W{z) 



SkHz)] 



Zi+Az/2 



-Az/2 



dz p{z - z'\az), 



(16) 



where S[k^-i,o,+i]{^) ~ [sinh(r), r, sin(r)], Az is the width 
of the bin, n{z) is the number density of source galaxies and 
7t,oo is the tangential shear for a source galaxy at z = oo, 
which is cancelled in taking the ratio of shears. 

p{z — z'\az) is the probability distribution for a galaxies 
position in redshift. If p{z — z'\az) is Gaussian then equation 
(16) simplifies to (Taylor et al., 2007; equation 25) 

= / dzn{z) " PAz[zi - z\az{zi)]Wiz). 

(17) 

For a sum of Gaussians each integral is simplified in the 
same way. The integrand PAz[zi — z\az{zi)] is given by (e.g., 
Ma et al., 2005) 



PAz[z\crz] 



erf 



erf 



z + Az/2 

z - Az/2 
V2a, 



(18) 



which is the part of the redshift error distribution that lies 
in a redshift bin centred on z of width Az, and erf (a:) is the 
error function. The function W{z) in equation (17) is given 

by 



W{Z) : 



w{z) 



dz'w{z')n{z')PAzlz - z'\crz{z)] 



(19) 



where 'w{z) is some arbitrary weighting function of the 
shears in redshift, which we take as 'w{z) = 1. 



Shesir-ratio intrinsic alignments 

Recasting the measured average tangential shear with an 
extra intrinsic shear component, 7i ^ ^i+Ci where d <C 7i 
the shear-ratio becomes 



7i_ ^ 7i + Si 



15 



7i 



I3 



(20) 



where the shear is the average in an aperture about a cluster 
in a given redshift bin. Since it is assumed that any intrinsic 
effect is small the shear-ratio signal is unaffected to first 
order by intrinsic alignment effects. The noise covariance 
properties of the shear ratio are affected in the following 
way. Taking the covariance of the shear ratio for two different 
pairs of bins (i, j) and (m, n) where Zi < Zj and Zm < Zn 
gives 



5Rij SR 

rr 



+ 



1 — 

7i7m 7j7n 

7»7" 7j7wi 
^77 

2,z,min(j,n) j-K 



fill 

2,max(i,m),3 j-K 



im,i>m 



CI 



n,j>n 



7j7n 



C2 



7i7rr 



7i7" 



/^7e /^7e 



+ 
+ 



7j7m Imlj 
y^ee flee 



7<7n 



7j7n 



7j7" 



(21) 



The first six terms are due to sample shear covariance, de- 
scribed in Taylor et al. (2007). The next four terms are the 
correlations between the intrinsic ellipticity of one bin with 
the tangential shear in another bin, correlations between ei- 
ther foreground bins or either background bins. Note that Cm 
is only correlated with 7, if Zi > Zm i.e. Cy^^^^^ = (dyidem) 
for all Zi > Zm- The next four terms, which have a negative 
constribution, are the correlations between the intrinsic el- 
lipticity of a foreground bin with the shear of a background 
bin CJ^ = (d'yjde-m)- The last four terms are the correla- 
tions between the intrinsic ellipticities between bins. Equa- 
tion (21) is as expected since there are no correlations be- 
tween the shear in a given bin and the intrinsic ellipticity in 
that bin i.e. {d'^idd) = 0. 

Since we use the tangential shear in an aperture about 
a cluster, equation (9) needs to be averaged. Assuming an 
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aperture of size 9 



2AGiE{zi,Zj) 

/ dee / de'e' 

Jo Jo 

f 

Jo 



422) 



Q2 + Q,2 _ 209' COS(0) + ( 

which can be calculated numerically, where E(zi,Zj) = 
D{zj)D{zi-Zj)/D{zi). We use values of Aqi = -0.92 x 10"® 
/iMpc~^ arcmin and 0o = 1-32 arcmin for tangential shear 
taken from Heymans et al. (2006). The II terms are given 
by 

All 

= l + (r[^,]-rN/S„)2 (23) 

where An = 0.45 x 10"=* and Bu = 1/t^^Mpc. Cf/ should 
be very small since Zi ^ Zj and the bin widths used, as a 
result of the photometric redshift error az{z) = 0.025(1 + ^:), 
correspond to Ar{z) = r[zj] — r[zi] ^ Bu. 

Note that since the intrinsic alignment effect does not 
enter into the signal of the geometric shear-ratio method 
the parameters Aqi and Oo cannot be marginalised over but 
rather a further noise term is added to the covariance. 



Shesir-ratio image distortion 

Since we assume a Singular Isothermal Sphere (SIS) profile 
(see Taylor et al., 2007) the shear can be written 



2\ey 



-i2<f> 



+ 7bK 



J, = _sR[7e-^*^^] = -^^7Cos(20)-7bu 



(24) 



where (pc is the azimuthal position of the a galaxy relative 
to the cluster center and (j) is the unknown systematic rota- 
tion of a galaxy's shape relative to the center of the galaxy's 
image. The shear-ratio signal takes the ratio of tangential 
shear so any overall distortion of the shear will cancel to 
first order, so that the signal remains unaffected. However 
the noise properties arc affected since the cross-component 
shear, which can be used to estimate the scatter in 74 (Kitch- 
ing et al., 2007) becomes 



(25) 



(7x)=(^)A7sin(20)+7bi, 



Using the notation of Taylor et al. (2007) the fractional error 
in the tangential shear due to the extended description of 
the shear is, to first order in 7biaB, 

2 j__.__2r 



A7» 
7i 



+ 



tan^(20) 
2iTnoAzi9^ 

tan^(2./>)7b ia: 
-KnoAzi9oA.y9 V,sin(2^) cos(2^) 



'KnoAziA^9'^ cos^(20) 
so that the overall fractional error in the shear is now 



^ 7i J 27rno Azi^o^^ cos(20) V 7» 



(26) 



(27) 



The effect of such the extra systematic terms is to in- 
crease the noise in the shear-ratio method. A caveat to this 
simple parameterisation is that any locally non-uniform dis- 
tortion, for example if A^ ^-f(^) oi 4> (t>{d), would not 
cancel in taking the ratio and any parameters describing 
such effects as a function of redshift or position would have 
to be marginalised over. 



APPENDIX B 

This Appendix will present in technical detail the effect of 
each primary systematic parameterisation on the 3D cosmic 
shear method. For a full exposition of the method see Castro 
et al. (2005), Heavens et al. (2006) and Kitching et al. (2007). 

We also take this oppurtunity to highlight an erratum 
present in Heavens et al. (2006) and Kitching et al. (2007). 
The signal covariance prefactor should not have a factor of 
Hq if k has units of Mpc~^. If k has units of /iMpc~^ then 
one should be present. This can be seen straightforwardly 
from Poisson's equation which enters the shear transform 
(equations 17 to 21; Heavens et al., 2006) 



2Ba 



S{k,l;r) 



(28) 



If k is in hMpc~^ then this should give no h dependence. 
Because Heavens et al. (2006) were using k in units of Mpc~^ 
this resulted in the transform being proportional to so 
that the covariance was proportional to h'^. 

Also in calculating the redshifts these papers used: Zi = 
Zi-i + {dz/dr)dr where dz/dr oc ft, but dr is oc h~^. So in 
calculating the Zi the factors of h should cancel. This was 
not done leading to an h dependence through the integrals 
J dz J dzp oc and a{r) = (1 + 2) oc (1 + ft) i.e. the signal 
covariance oc ft®. So that the overall spurious ft dependence 
of the signal covariance was a factor of ft^'' overall. 

These errata resulted in the constraints on ft being too 
optimistic in these papers through this strong ft dependence. 
Fortunately, for the fiducial surveys considered in these pa- 
pers, the ft constraint is dominated by the CMB Planck 
prior, so that there is a small effect on the dark energy par 
rameter constraints. For the larger surveys considered in this 
paper the spurious ft dependence would lead to the lensing 
constraint becoming better than the CMB Planck prior and 
as such would effect the dark energy constraints by lifting 
parameter degeneracies. 

3D cosmic shear photometric redshift systematics 

3D cosmic shear uses the 3D spherical harmonic coefficients 
characterised by azimuthal £ and radial fe-modes {k is in 
ftMpc"'^) given by 



%{k,£) = '^■yikje{kr^)e 



(29) 



for a galaxy at position (r^, 9^) with shear 7^, ji{z") are 
spherical Bessel functions. Note that this results in four in- 
dependent data vectors ^i{k,£), ^l{k,£), ^2{k,£), ^i{k,£) 
where -R denotes a real vector and / the imaginary part. 

Since the mean of such coefficients is zero it is the covari- 
ance of the coefficients that is used to extract cosmological 
information. This covariance is a sum of signal and noise 
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terms C = S + N where the signal is given by (Heavens et 
al., 2006; equations 23 to 26) 

s = (f{k,£)f*ik',e'))s = Qi{k,k')s''ie-e') (30) 

where Qe {k, k') can be written as 



■J ^Ge{k,k)Gi{k',~k) (31) 



Qe{k,k') 
where 

Ge{k,k) = k J dzdZpn^{zp)W{zp)p{zp\z)Ue{r,k)je{kr'') 

(32) 

where p{Zp\z) is the photometric redshift probabihty distri- 
bution and 

Ui{r, k) = r df ^"^^:^^ ^Ps{k;f)n{kf) (33) 
Jo '^V) 

where r = r{z) and 

FK{r,r') = {Su{T-T')/[Sk{r)Sk{r')]] 



(34) 



Qrn is the present dimensionlcss matter density, Hq is the 
present Hubble parameter, P5{k; r) is the matter power spec- 
trum, n{z) is the source number density distribution and 
W{z) is an axbitrary weighting function which we set to 
W{z) = 1. 

The photometric redshift probability distribution enters 
the covariancc in equation (32). It acts to damp the signal ra- 
dially at scales approximately equal to and loss than the pho- 
tometric redshift error i.e. at k values of fc > 27r/i/(3000a2). 
Adding extra systematic parameters to p{zp\z) adds extra 
parameters to the signal of the 3D cosmic shear method but 
no extra covariance. 



3D cosmic shear intrinsic alignments 

Intrinsic alignment effects are included in the 3D cos- 
mic shear method by including additional signal covariance 
terms C={S + GI + II) + N. 

The effect of the GI intrinsic alignment effect for 3D 
cosmic shear is to add a further covariance between shear 
and intrinsic ellipticity (■y{k,£)e*{k' ,£')). Using the notation 
of Heavens et al. (2006) , including a convolution over pho- 
tometric redshift, this can be written as 

GI = {^{k,£)e*{k\l')) 

= ^^~3^ j dzdZpd^O j dz'dz'pd^d' 
p{z\zp)p{z'\z'p)n{z)n{z')kk' 

k{kr\z])3t{k'r\z'])e-'^-^e+'^' {^{r)e*{r')). (35) 

Substituting the parameterised form from equation (9) this 

can be written as 

{^{k,£)e*{k',£')) = [0^ j dzdz' 

E{z,z')Hi{z,k)He{z',k')I{n) (36) 
where E{z,z') = D{z')D{z - z')/D{z) is the lensing effi- 



ciency function and 

He{z,k) = J dzpp{z\zp)n{z)je{kr[z])k 



(37) 



m) = .^..^.'e-<^-^'-)^^^(38) 

where Q. = AO x is the solid angle of the survey. The 
integrals in /($!) can be solved numerically. Wc use values 
of ^Gi = — 1.26x 10~^ /iMpc~^ arcmin and = 0.90 arcmin 
taken from Heymans et al. (2006). 

The II intrinsic effect is included by using the following 
covariance, including a convolution over photometric red- 
shift 

II = {e{k,e)e*{k',£')) 

= (i^) / dzdz'M{k,z)M{k',z') 
1 



^ rlz]-rlz'] y 



(39) 



where O is the solid angle of the survey, r{z) are comoving 
distance and 

M{k,z)= J dzpp{z\zp)n{z)je{kr)k. (40) 

We use values of An = 0.45 x 10~^ (for mixed galaxies) and 
Bu = lft~^Mpc taken from Heymans et al. (2006) 

The intrinsic alignments systematics add further covari- 
ances to the 3D cosmic shear signal and the extra parameters 
An, Agi and 9o- 



3D cosmic shear distortion 

Transformed shear components (71,72) can be written using 
the parametcrisation of equation (8) as 



7^ = ^7 [71 cos{2(j)) + 72 sin(20)] + 7ibias 
72 = yl-^[-7i sin(20) -|-72COs(2c/>)] 4- 72bia 



(41) 



where 7bias = 7ibias + *72bias- So that, for example, the new 
72 estimator ^2{k,t) can be written, using the notation of 
Heavens et al. (2006), 

1/2 

l2{k,£) = 



E 



k{A^[-^-i sin(20) -|- 72 cos(20)] -|- 72bias} 



ji[kr{z)]e 



il.e 



(42) 



Taking the covariance results in four terms respectively pro- 
portional to 

4(7i7i*>sin2(2<^) 



^7(7272) cos^(2(/>) 

^7(72bias72) COs(2(?i) 
A(7272biaB)cOs(2</>) 



(43) 



to first order in 72bias- 

In the calculation, in order in save computational 
time, the 72*(fc,£) component is chosen as representative 
so that in the signal covariance, equation (31), Xj^ = 
2\£\^ cos{(f>e) sm{(j)e) where (pe is the angle of the £ vector in 
the {£x, iy) phase space. We only consider > to avoid 
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double counting. By choosing (pi = 7r/4 the prefactor be- 
comes X£ = the Fisher matrix is then integrated over 
all modes in a given shell using the 7^(fc,^) component as 
representative 

Fci3=g [ ' d4>i f dUF^nil) (44) 

J-it/2 J 

where F„/3(£) = \T!t\^{C^)-''C%{C^)-''C^^. The density 

of states in ^-space due to the survey size is g = Af2/(27r)^ 
so that 

Fc0 = ^ J d££Fc0{i). (45) 

To test this approximation a full Fisher matrix calculation 
was done over all (j>e where Fapit) —> Fapil-, 4>i), for all four 
data vectors, the different (f>e dependence comes from the 
prefactor. The dark energy parameter errors were found be 
in agreement to within ±0.001 and since the computational 
time is 4 X A^^ larger for the full calculation (where A^^ is the 
total number of modes calculated in ^-space) the approxi- 
mation is used for all predictions in this paper. 

Since in the calculation we use 0^ = 7r/4 as a represen- 
tative point in the £ phase space and integrate over £ the 
terms in equation (43) can be reduced since (7171) 0. 
Also, since we assume that 72bias =constant the expres- 
sion reduces further since (72bias72) = 72bias(72) = 0. 
So the signal is simply multiplied by an extra factor of 
Aj cos^(2(^). In the shot noise calculation {7171) — {7272) = 
since [{7171*) sin2(20) + (7272) cos2(2(/))] = [cos^{2(l>) + 
sin^(2(;6)](7e = iTe, so the shot noise is simply multiplied by 
A^. A caveat is that the assumption that the mean signal 
is zero would break down if 7bias became too large. Here we 
assume that the assumption that the mean signal is zero is 
still valid since 7bias <S 7 which will be valid to first order 
in 7bias- 

Since the extra shear systematic parameters are in 

the signal, the calibration A-., and the angle d> can be 
marginalised over and no extra covariance terms are added 
to the method. 
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